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DIFERENCIA CION NUMERICA 


La derivada de una función tiene muchas aplicaciones, entre las cuáles esta la 
determinación de la velocidad instantánea de una partícula o móvil a partir de su 
función de posición. Este proceso es en ocasiones algo muy sencillo cuando se cuenta 
con dicha función, pero cuando se requiere solucionar el mismo problema con un 
conjunto de datos discretos y no con su función, el procedimiento no puede ser llevado 
de igual manera, es decir, el calculo no nos da una solución directa, por lo tanto se debe 
recurrir a otro tipo de análisis. 

En la unidad anterior se estudio la aproximación polinomial para un conjunto de datos 
tabulados con base en esos conceptos. En esta unidad se estudiara la aproximación 
polinomial a partir de algunos valores de una función f, y dando ciertos puntos se puede 
determinar su derivada en un punto dado. 

La derivada de una función/ enxo, esta definido de la siguiente manera: 


/'(x 0 ) = lim 

/ î —>0 


/(x 0 + h) — /(x 0 ) 
h 


Para valores pequenos de h, podemos aproximar la derivada de / en xo de la siguiente 
manera: 


/’(x 0 ) 


/(X o +h)-f(x 0 ) h ^ Q 
h 


(4.1) 


Notaremos la aproximación a la derivada de una íunción/corno/’ 

Esta formula aunque sencilla no tiene un comportamiento estable, ya que para funciones 
lineales puede llegar a ser exacta, no siendo así para funciones más generales. Pero sin 
duda alguna, es un buen punto de partida para el calculo de la derivada de una función, 
además hay que considerar que en algunos casos es la única opción con que se cuenta. 

Para estimar el valor del error asociado a la ecuación (4.1), nos valemos del polinomio e 
Taylor de primer grado cuya estructura es la siguiente: 
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/O) - f(x o) + /'O 0 )(x - x 0 ) + (X - x 0 )■ 2 (4.2) 

Si se toma x = x 0 + h ,entonces h = x - x 0 y se reemplaza en el polinomio (4.2), se tiene 
entonces que: 


/(x 0 +h) = f(x 0 ) + /' (x 0 )h + 



x 0 < Ç < x 0 + h 


Despejando f’(xO), se obtiene 


/'(, 0 ) = ++A'+/+>_ o; , 

h 


(4.3) 


f"(â) 

Donde O = - -—, es el error de truncamiento. 

2 

Si h> 0, a la formula (4.3) se le denomina la primera diferenciafìnita hacia delante o 
diferencia progresiva También podemos obtener la diferencia finita hacia atrás o 
diferencia regresiva si h < 0 


/'Oo) 


/Oo)-/Q 0 -h) { oh 

h 


(4.4) 


Ahora bien, si sumamos las ecuaciones (4.3) y (4.4), obtendremos la diferencia fìnita 
centrada de la siguiente manera: 


0 2 h 


(4.5) 


Podemos observar que el error de la ecuación (4.5) es del orden h 2 , a diferencia de las 
ecuaciones (4.3) y (4.4) que tienen un error del orden h, es decir que la ecuación (4.5) 
converge rápidamente a cero, pero para ello se debe contar con 3 valores de f(x) a 
diferencia de (4.3) y (4.4) que solo requiere de dos puntos de la función f(x). 
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Fig. 4.2 Diferencia Finita Regresiva 
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De manera análoga a la interpolación polinonrial, el uso de mas puntos en la evaluación 
de la derivada producirá mayor exactitud; aunque esto inrplica mayor cantidad de 
evaluaciones funcionales y aumento de error de redondeo. 


Entre las formulas mas conrunes están la de tres puntos y la de cinco puntos. A 
continuación se nrostrara la deducción de la fornrula de tres puntos a partir del 
polinonrio de Taylor de segundo grado. 

Se sabe que: 


f(x) = f(x 0 ) + f(x 0 )h + 


f"(x 0 ),2 


/"'(£> 


Despejando f(xo) se tiene 


f'(x 0 ) 


f(x)-f(x 0 ) f"(x 0 ) h f"(Ç ) h 2 

h 2 6 


Conro se quiere detemrinar la fornrula de los tres puntos, se tonra 
XO, 

Xi = Xo+ h , 


X 2 = xo+2h 

y usando la ecuación (4.6) x = xi 
f(x 0 +h)~f(x o) 


f'(x 0 ) = 


h 


r+A h J+Sïe 
2 6 


(4.6) 


Para determinar f”(xo) en térnrinos de xO, xl y x2, se realiza los polinonrios de Taylor 
de segundo grado para X/ =x 0 + h y para x 2 = xo+2h 

f(x 0 + h) = f(x 0 ) + f'(x 0 )h + f " ( 2 X(>) h 2 (4.7) 

f(x 0 + 2/0 = f(x 0 ) + / '(x 0 )2 h + í ^l(2h) 2 (4.8) 

nrultiplicando a (4.7) por -2 y sumándoselo a (4.8) se tiene: 

f(x 0 + 2/0 = f(x 0 ) + f \x 0 )2h + 2f"(x 0 )h 2 
- 2 f(x 0 + h) = -2f(x 0 ) - / \x 0 )2 h + f"(x Q )h 2 


f(x 0 + 2/0 - 2 f(x 0 + h) = -f(x 0 ) + f"(x 0 )h 2 

Despejando/’ ’(x 0 j de la ecuación (4.9) se tiene: _ 
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= /(J °~ 2/fa +/»+/(*.+2 '») +0h 
h 

Reemplazando (4.10) en (4.7) se obtiene 

fl(x , = f(x 0 +h)-f(x 0 ) f(x 0 )~ 2f(x 0 + h) + f(x 0 + 2 h) h + ^ 

J 1 0 h 2h 2 

r(x _ 2f(x 0 + h)~ 2 f(x 0 ) - f{x Q) + 2/(x„ +/Q- /(x 0 + 2h) ^ Qh 2 


/■fa) = -Vfa) + 4/fa +ft)-/fa +2<.) +0;<î (4U) 

2 h 

La ecuación (4.11) es conocida como la formula de los tres puntos progresiva. De 
manera similar es posible desarrollar las formulas de los tres puntos regresivas y 
centradas, así corno las de los cinco puntos y derivadas de orden superior. 

Las tablas a continuación resumen las formulas para aproximación de las derivadas 


Primera Derivada 


Formula de 

diferencias 

finitas 

progresivas 

Dos Puntos: 

r(x j_f(xo+h)-f(x 0 ) 
h 

Oh 

Tres Puntos: 

r(x , _ _ 3/(x 0 ) + 4/(x 0 +h)-f(x 0 +2h) 

J 1 0 2 h 

Oh 2 

Formula de 

diferencias 

finitas 

centradas 

fu x \ _ f( x o +h )-f( x o~h) 

J K 0 2 h 

Oh 2 

, f(x 0 -2h)-%f(x 0 -h) + %f(x 0 +h)-f(x 0 +2h) 

Oh 4 

/(Xo)_ 12 h ^ ^ 

Formula de 

diferencias 

finitas 

regresivas 

j; (x ) _f( X o)-f( X 0- h ) 
h 

Oh 

fu x ) _ f ( x 0 -2/0 -4 f(x 0 -h) + 3 f(x 0 ) 

J 1 0 2 h 

Oh 2 


Segunda Derivada 
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Formula de 

diferencias 

finitas 

progresivas 

f ,,,„ , _ f(x 0 )~ 2 f(x 0 +h) + f(x 0 +2h) 

J V x 0 / — J 2 

n 

Oh 

2 f( x o)~ 5/(x 0 + h) + 4 f(x 0 + 2h) - f(x 0 + 3 h) 

J “ j 2 

n 

Oh 2 

Formula de 

diferencias 

finitas 

centradas 

,,,,„ , _ f(x o - h)- 2/ (x 0 ) + f(x 0 + h) 

J V*0/ — J 2 

n 

Oh 2 

f ,,,„ \ _-/(*b”2A)+16 /(x 0 -/ 2 )- 30 /(x 0 )+ 16 /(x 0 +h)-f(x 0 +2 h) 

Olì 4 

' ° \2lr 

Formula de 

diferencias 

finitas 

regresivas 

S,,,„ \ _ f ( x 0 - 2 / 2 ) - 2/(x 0 -h) + f(x 0 ) 

J V x 0/ — J 2 

n 

Oh 

f»í~\_-f ( x o - 3/0 + 4/ (x 0 - 2 h) - 5/ (x 0 - h) + 2/ (x 0 ) 

J V X 0/ J 2 

n 

Oh 2 


Tercera Derivada 


Formula de 

diferencias 

finitas 

progresivas 

~ f( x o) + 2 f( x o + h)~3f(x 0 +2h) + f(x 0 + 3 / 2 ) 

J v-^o / — J 3 

n 

Oh 

„,,,„ \ _ -5/(a) +1 Wo +h)~ 2 4f(x 0 +2h)+l4f(x 0 +3h)-3f(x 0 +4h) 

J (V 2^3 

Oh 2 

Formula de 

diferencias 

finitas 

centradas 

f ( x 0 - 2/0 + 2 f(x 0 - h) - 2 f(x 0 + / 2 ) + f(x 0 + 2 / 2 ) 

1X0 1 8/2 3 

Oh 2 

\ fh -30-8/Cy, -20+1 J(s> ~h)~\ f(s } +h)+%f(Xo +2h)-f(s ) +3h) 

Oh 4 

' W 8 ìf 

Formula de 

diferencias 

finitas 

regresivas 

f,,,,„ \_ ~ f ( x 0 - 3 / 2 ) + 3/(x 0 - 2 / 2 )-3/(x 0 -h) + f(x 0 ) 

J v-^o / — 7 3 

n 

Oh 

f ,„, \ 3/0)-4/2)-14/tc-3/2)+24/(xq-2/;)-l §f(^-/2)+5/(0 
' W 21ì 

Oh 2 


Ejemplo 4.1 

Aproximar la primera derivada de la función f(x) = e 2x con todas las formulas que se 
encuentran en las tablas para unxs= 1.1 con un tamano de paso h = 0.1. Use el hecho 
que f’(l.l) = 18.050 para determinar el error absoluto y relativo producido por cada 
formula. 
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Solución: 

Usando las formulas de diferencias finitas progresivas (tabla 4.1), se requieren los 

valores de f(xo), f(xo+h) y f(xo+2h) 

f(x 0 ) : f(U) = 9.025 

f(x 0 +h) : f(1.2) = 11.023 

f(x 0 +2h): f(1.3) = 13.464 


/'( 1 . 1 ) 


11.023-9.025 

0.1 


19.980 


E a = |18.050-19.980|= 1.930 
E r = 1.930 / 18.050 x 100% = 
10.7% 


/'d.l) 


-3-(9.025) +4-(11.023)-13.464 ?65 E A = 118.050-17.765| = 0.285 

2 (°.l) E r = 0.285 / 18.050 x 100% = 

1.58% 


Usando las formulas de diferencias finitas centradas se requiere: 


f(x 0 -2h): f(0.9) = 6.050 
f(x 0 -h) : f(1.0) = 7.389 
f(x 0 +h) : f(l.2) = 11.023 
f(x 0 +2h): f(1.3) = 13.464 


/'(%) 


11.023-7.389 

2 ( 01 ) 


18.170 


E a = |18.050-18.170| = 0.120 
E r = 0.120 / 18.050 x 100% = 
0.665% 


/'(%) 


6.05 0-8-(7.3 89)+8-(11.023) -13.464 
12 ( 0 . 1 ) 


= 18.048 


E A = |18.050-18.048| = 0.002 
E r = 0.002 / 18.050 x 100% = 
0.009% 


Ahora usando las diferencias finitas regresivas, se requieren los siguientes valores 


Ing. Jorge Mario Pena Consuegra 
Ing. Marvin Enrique Molina Cárdenas 
7 




MÉTCmOS nvmbrjcos 


f(x 0 -2h): fïfì.9) = 6.050 
f(xo-h) : f(1.0) = 7.389 
f(xo) : f(U) = 9.025 

E a = |18.050-16.360|= 1.69 
E r = 1.69 / 18.050 x 100% = 
9.363% 


/'(%) = 


9.025-7.389 

0.1 


= 16.360 


/'(%) 


6.050 -4(7.389) + 3(9.025) 

2 ( 0 . 1 ) 


17.845 


E a = |18.050-17.845| = 0.205 
E r = 0.002 / 18.050 x 100% = 
1.136% 


Del ejercicio anterior se puede concluir que las mejores aproximaciones las ofrecen las 
diferencias finitas centradas, por lo tanto estas se deben preferir siempre que se puedan 
utilizar. 
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DIFERENCIACIÓN NUMÉRICA CON MATLAB 

Para el calculo de la derivada de una función a partir de valores tabulados de la función 
se realizo el siguiente archivo .m. 

Este archivo .m encuentra la derivada para un conjunto de puntos igualmente espaciados 
en x. Estos puntos se deben introducir en forma de vectores. Un vector para x y un 
vector para y. E1 archivo .m genera otro vector d que es un vector que contiene la 
derivada de la función dada con x e y encontrado usando la formula de los dos puntos. 
Para los puntos extremos el programa utiliza los métodos de diferencia finita progresiva 
(para el primer punto) y regresiva (para el primer punto). Para los puntos centrales se 
utiliza el método de diferencia finita centrado 


der2p.m 

function d=der2p(x,y); 

%d=der2p(x,y) 

%archivo .m que sirve para determinar la derìvada de una función dada 
por 

%medio de puntos igualmente espaciados en x con tamano de paso h, por 
medio 

%de la formula de dos puntos. 

%En los puntos extremos se utilizan las formulas de diferencias 
fìnitas 

%progresivas y regresivas y en los puntos medios se usa la formula de 
%diferencias finitas centradas. 

%las variables x y y son vectores de igual tamano que ìncluyen los 
valores 

%de la función en x y en y. Estos vectores deben ser de tamano mayor 
o 

%igual que 3 


ìf size(x)==size(y) %se observa que los dos vectores sean del mismo 
tamano 

[m, n]=size(x) ; 

h=x(2)-x(l); %tamano del paso 

if n>=3 %se verifica que se den mas de dos puntos 

for i=2:n-l %se trabajan con los puntos ìntermedios 
d (i) = (y (i + 1) -y (i-1) ) / (2*h) ; 

end 

d(1)=(y(2)-y(1))/h; 
d(n)=(y(n)-y(n-l))/h; 

else 

d='se deben dar mas de dos puntos' 

end 

else 

d='los vectores x y y deben ser del mismo tamano' 

end 


Ing. Jorge Mario Pena Consuegra 
Ing. Marvin Enrique Molina Cárdenas 
9 





métovos mJMmjcos 


Encontremos la derivada de la función f(x) = e 2x desde x = 0.8 hasta x = 1.3 con un 
tamano de paso h = 0.1, usando este archivo .m 


COMMAND WINDOW 

EXPLICACION 

» x=0.8:0.1:1.3; 

Se crea un vector con elementos 

linealmente espaciados que comienzan en 

0.8 y termine en 1.3 con un tamano de paso 

h = 0.1. E1 punto y coma es para que no se 

visualice el vector x. E1 vector x es de 6 

elementos. Si quiere verlo escriba x y 

oprima enter 

» y=exp(2*x); 

Se evalúa la función f(x) = e Zx , usando 

para ello el vector x. Se recibe un vector 

de igual tamano que el vector x donde cada 

elemento de y es la evaluación de la 

función usando cada elemento del vector x 

» d=der2p(x,y) 

Se utiliza el archivo der2p.nr para 

determinar la derivada en cada uno de los 

puntos dados. Se crea el vector d que se 

muestra a continuación. 

d = 

10.9662 12.1801 14.8768 18.1706 22.1936 24.4056 
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Se obtiene entonces la siguiente tabla de datos: 


X 

y 

m 

0.8 

4.953 

10.966 

0.9 

6.0496 

12.18 

1 

7.3891 

14.877 

1.1 

9.025 

18.171 

1.2 

11.023 

22.194 

1.3 

13.464 

24.406 
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